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Towards a non-perturbative matching of HQET and QCD with dynamical light quarks 



1. Introduction 

In the light of the expected progress in flavour physics thanks to the impending B -physics 
experiments [|l], 0, precision lattice QCD more and more becomes to play a crucial role for a quan- 
titative and accurate interpretation of these experimental results in the framework of the Standard 
Model and beyond, since it provides a theoretically sound approach to non-perturbatively compute 
the contributing matrix elements of operators among hadronic states. 

A particular problem of dealing with heavy-Ught systems involving the b-quark as the heavy 
flavour by means of lattice QCD consists in the two disparate intrinsic scales that actually accom- 
pany any lattice calculation: the lattice spacing, a, has to be much smaller than 1/nib in order to 
allow for a fine enough resolution of the B -meson states in question, and the linear extent of the 
lattice volume, L, has to be large enough for finite-size effects to be under control. Heavy Quark 
Effective Theory (HQET) on the lattice however, which relies upon a systematic expansion 

of the QCD action and correlation functions in inverse powers of the heavy quark mass (m) around 
the static limit, offers a formally reliable solution to this problem. Still, for lattice HQET and its 
numerical applications to lead to precise results with controlled systematic errors in practice, two 
shortcomings had to be left behind first. 

One is the exponential growth of the noise-to-signal ratio in static-light correlation functions, 
which is a consequence of the appearance of power divergences in the effective theory. As demon- 
strated in studies in the quenched approximation [|5]-^ as well as in the theory with Nf = 2 dy- 
namical quarks JIol], this problem can be overcome by a clever modification of the traditional 
Eichten-Hill discretization of the static action. 

Another difficulty, more serious on the theoretical level, is associated with the aforementioned 
power divergences. Since in the effective theory mixings among operators of different dimensions 
are present, already the static limit of HQET is affected by a power-law divergent (~ g^/a) additive 



mass renormalization. Unless the theory is renormalized non-perturbatively []11|], it follows from 
this power-law divergence of lowest-order HQET — and, of course, from further ones ~ g^/a"'^^ 
that arise at 0(l/m"), « > 1 — that the continuum limit does not exist owing to a remainder, which. 



at any finite order [|12|-|14p in perturbation theory, diverges in the continuum limit. 

In ref . [ 15 1 a general solution to the latter has been worked out and numerically implemented 
for a determination of the b-quark's mass in the static and quenched approximations as a test case. 
The method is based on a non-perturbative matching of HQET and QCD infinite volume. It was 
subsequently extended to also include the 0(l/m) terms into the quenched computations of the 
mass of the b-quark, m^^(mb) JT^ (see refs. [T^ for recent reviews in broader context), and 
of the B-meson decay constant [|I9|]. 

An attractive property of the strategy, briefly summarized in Section ||, is that most parts of the 
actual calculation do not involve very large lattices. Hence, it is natural to remove the quenched 
approximation as the dominating remaining systematic uncertainty in our previous works using 
this method. The additional computational effort required if dynamical quarks are included is 
only moderate, except for the last step that involves the extraction of B-meson properties from 
simulations in physically large volumes (with spatial extents of ?a 2 fm or more) and thus will be 
computationally much more demanding than the finite- volume simulations for the non-perturbative 
renormahzation part. 
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In the present report we outline the various steps towards an HQET computation of the mass of 
the b-quark including the 0(1 /m) correction along the lines of refs. [|^, [T6|] in two-flavour lattice 
QCD, where most of the emphasis is put on the renormalization of the effective theory through 
the non-perturbative matching to QCD in finite volume in order to perform the power-divergent 
subtractions. This step requires, in particular, a determination of the relation between the renormal- 
ization group invariant (RGI) and the subtracted bare heavy quark mass in the relevant parameter 
region of A'f = 2 QCD, which we present together with numerical results on the corresponding 
renormaUzation constant and improvement coefficients in some detail in Section ^. Results from 
the matching itself, which has just been started at the time of writing, as well as from the necessary 
simulations of the effective theory in small and intermediate volumes will only be available at later 
stages of our project. 

2. Survey of the computational steps 



HQET 
in large 
volume: 



TTLb 



0.5 fm 



matching 



0"k 



Figure 1: The strategy for performing computations in lattice HQET via a non-perturbative determination 
of the HQET parameters from QCD simulations in a small volume. It is designed such that steps indicated 
by arrows are to be repeated at smaller lattice spacings to reach a continuum limit. 



Let us briefly recall the general strategy, introduced in []15|]. It allows for a formulation of 
(zero- velocity) HQET in the framework of lattice QCD, where all steps of the computations includ- 
ing the renormalization are carried out non-perturbatively and the continuum limit can be taken. 

The basic idea is illustrated in figure [I] and starts from a finite volume of extent L\ ^ 0.5 fm. 
There, one chooses lattice spacings a sufficiently smaller than 1/mb such that the b-quark propa- 
gates correctly up to controllable discretization errors of order a^. Since the relation between the 
RGI mass and the bare mass in QCD is known pO|], suitable finite- volume observables ^k{Li,M\,) 
can be calculated as a function of the RGI b-quark mass, Mb, and extrapolated to the continuum 
limit. The next step is to perform the power-divergent subtractions non-perturbatively by a set of 
matching conditions, in which the results obtained for are equated to their representation in 
HQET (r.h.s. of figure |]). At the same physical value of Li but for resolutions L\/a = 0(10), the 
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previously computed heavy-quark mass dependence of <I>i.(Li,Mb) in finite-volume QCD may be 
exploited to determine the bare parameters of the effective theory for a k, (0.025 — 0.05)fm. In 
order to evolve the HQET observables to large volumes, where contact with experiments can be 
made, one also computes them at these lattice spacings in a larger volume, L2 = 2L\. The resulting 
relation between <I>j-(Li) and (^^{1-2) is encoded in associated step scaling functions Ok, as indicated 
in figure D as well. Finally, the knowledge of <I><:(L2,Mb) and employing resolutions L2/a = 0(10) 
fixes the bare parameters of the effective theory for a ^ (0.05 — 0.1)fm so that a connection to 
those lattice spacings is established, where large-volume observables, such as the B-meson mass 
or decay constant, can be calculated (l.h.s. of figure |T]). 

Having in mind the computation of M\, as the specific application, this sequence of steps yields 
an expression of wib (taken to be the physical input) as a function of Mb via the quark mass depen- 
dence of ^k{Li,M\,), which eventually can be inverted to arrive at the desired physical value of 
the RGI b-quark mass extracted from the effective theory. As pointed out before, the whole con- 
struction is such that its various pieces separately have a continuum limit. With the realization 



of this strategy for the quenched case it was shown in ref. [ |16| ] that a determination of Mb in- 
cluding 0(l/m) in HQET only requires up to three matching observables, Oi, <I>2 and O3, if the 
spin-averaged B-meson mass is used as physical input. That is also the path we will follow in our 
present extension to the case of two-flavour QCD. 

2.1 Definition of the matcliing volume 

We consider QCD with Nf = 2 mass-degenerate dynamical quarks, which are identified with up 
and down. All other quarks are treated in the quenched approximation. A particularly convenient 
renormalization scheme, in which finite- volume observables suitable for a non-perturbative match- 
ing of the effective theory with QCD can readily be constructed [ [l5l ^Aj 22], is the Schrodinger 



functional (SF) []23|]. Relativistic and static quarks were introduced in [ ]24| ] and []25|], respectively. 



where in the latter reference it was found that the HQET expansion of the boundary quark fields is 



trivial up to and including 1 /m-terms.^ Adopting any unexplained notation from refs. [15, 25], we 
only mention the periodicity phase 6 of the fermion fields as a further kinematic parameter and the 
fact that homogeneous Dirichlet boundary conditions in time at ;co = and xq = T are employed. 
Since the parameters 6 and masses of the quenched quarks can be set independently of those of 
the sea quarks, the basic situation for extracting heavy-light physics from SF correlation functions 
is the same as in the quenched approximation 1 15, |l^]. Moreover, in the finite-volume simulations 
we set 6 = 0.5 for the dynamical light quarks and their PCAC mass to zero, mi = 0. 

The quantities that enter the non-perturbative matching procedure described above have 
to be evaluated in the continuum limit. To this end we want to compute them for a series of 
bare parameters {L/a,l5,K\) such that the renormalized parameters in the light quark sector are 
fixed and thereby physics is kept constant along the approach to the continuum limit. Here, Kj 
denotes the hopping parameter of the dynamical light quarks. Our constant physics condition on 
the renormalized SF coupling, g^{L), and the light quark mass reads 

g2(Lo) =2.989, ^o = y, Lomi(Lo)=0. (2.1) 



'From now on, m generically denotes the mass parameter of the heavy quark treated in the effective theory, while 
the masses of the non-degenerate quark flavours in the relativistic theory are distinguished explicitly where necessary. 
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This choice now defines the spatial extent Li of the volume, in which the matching between HQET 
and QCD is performed. Although an exact knowledge of Li in physical units in not yet needed at 



this stage, one can already infer from the known running of the SF coupling for Nf = 2 [ [26| | that 
L\ 0.5 fm. Hence, we will finally have L2 = 2L\ f» 1 fm and thus Loo = 4Li f« 2fm for the large 
volume, which is well consistent with the envisaged strategy, figure |l|. 

We have fixed f^(Li/2) = 2.989 by a new simulation at Lq/a = 20, T = Lq, and made tentative 
interpolations in j8 = 6/go for given LQ/a< 16 to this target value, based on the known dependence 
of the SF coupling and the current quark mass on the bare parameters (j8 , k) available from the data 
of ref. [p6|]. Using the known j8-function and our experience from the quenched calculation []T^], 
we can estimate that an uncertainty of about 0.04 in the coupling will translate via the resulting one 
in Li into an uncertainty in the b-quark mass of at most 0.5%. The condition of zero light quark 



mass in eq. (gJJ) is met by setting k: = Kj to the critical hopping parameter, jq, estimated again 



on basis of published data []26|], whereby a shght mismatch of \LQm\{LQ) \ < 0.05 of this condition 



is tolerable in practice. The triples {Lo/a,[5, Kj), which approximately define the extent Li of the 



matching volume through eq. (|2JJ) and which are used in our subsequent study of improvement 
and renormalization factors, are collected in columns 2 - 4 of table [l| in Section ^. 

The preliminary interpolation procedure for (j3 ) underlying these j8-values is currently be- 
ing checked (and refined) by direct simulations, in order to avoid a non-negligible systematic error 



from small violations of the condition (2.1) on the final results. Yet, this will affect our estimates 
of bm and Z in Section ^ only at a negligible level, because there any deviation from the line of 
constant physics only entails a small change of the O(a^) effects. 

2.2 Fixing the heavy quark mass in finite-volume QCD 

Having fixed Li via enforcing constant physics at Lq = L\/2, the computation of the heavy 
quark mass dependence of the finite- volume observables (^k, which is the key element in the non- 
perturbative matching step within our strategy, will amount to evaluate heavy-light SF correlation 
functions in a volume L\xT , T = L\, for a series of precisely fixed values of the renormalized 
heavy quark mass covering the b-quark mass region. 

This is achieved by exploiting the O(a') improved relation between the (subtracted) bare heavy 
quark mass mq h and the RGI mass [ pO[ , viz. 

M = h{LQ)Zm{go,LQ/a) m^^h (l +^m(,go)a'«q,h) + O (a^) , (2.2) 

where 

^' ?0)ZA{g0) If I 1 

-— , arria h = ~ 

Zp{go,Lo/a) ^' 2\Kh Kc 

and Za is known non-perturbatively from ref. p8|]. The scale dependent renormalization constant 

Zp may be calculated for the relevant couplings on Lq x T lattices with T = Lq in the same way as 

in ref. [pO|]. The factor 

MZ.o) = =^ = 1.521 (14), ^ = L = ^, (2.4) 

represents the universal, regularization independent ratio of the RGI heavy quark mass, M, to the 
running quark mass in in the SF scheme at the renormalization scale /Xq. h{Lo) was evaluated by a 
reanalysis of the Nf = 2 non-perturbative quark mass renormalization data published in ref. [pO|]. 
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Therefore, in order to specify the hopping parameters of the heavy flavour, Kj,, which according 
to eq. ( 12 ) for given L\/a = 2Lq / a and j8 = 6/gQ correspond to a series of dimensionless RGI quark 



mass values z = LiM in the b-quark region, it remains to accurately determine the improvement 
coefficient and the renormalization constant Z. We discuss this computation for the relevant 
weak coupling range (cf. table [I]) of 0(a) improved two-flavour lattice QCD in the next section. 

2.3 Preparing for the finite-volume computations 
2.3.1 Matching to QCD 

As mentioned in the foregoing subsection, on the QCD side this step consists in calculating the 
quark mass dependence of the quantities <I>i, <t>2 and <I>3 in the volume Lj. For the exact definitions 



of these effective heavy-light meson energies in terms of SF correlators we refer to ref. []16|]. In 
addition to Lj x 7, J = Li, lattices with T = L\/2 will also be needed (cf. Appendix C of [|l6|]). 

The aforesaid fine-tuning of j8 for L{/{2a) < 16 to satisfy the condition g'^{L\/2) = 2.989, 
eq. (2.1), with a precision Ag^ < 0.04 requires up to (Li/a)^ = (2Lo/a)'* = 40^* lattices with A/^f = 2 



at sea quark parameters close to those quoted in table [T] in order to reach the continuum limit. 

2.3.2 Parameters for HQET simulations inL^ x T with T = Ly,Li/2 

For the determination of the step scaling functions Li —>■ L2 = 2L\ belonging to the <I>i:'s coun- 
terparts in HQET, we must fix the simulation parameters for resolutions 6 <Li/a < 16. The cor- 
responding constraint on the renormalized coupling at Li is g^{Li) = a(2.989) = 4.484(48) [^6|]. 

As a starting point for the tuning of g^{Li) at each Li/a, we introduce another low-energy 
scale, L*, defined via g^(L*) = 5.5 and obeying [ |29| ] 



In (LVa) = 2.3338 + 1.4025 (j8- 5.5) , jS G [5.3 , 5.8] , L7a G [7.8 , 16.1 ] , (2.5) 

which allows to estimate the ratio ri = L\/L* ^ 0.8 in the continuum limit. Trial j8-estimates for 
the range of Li/a in question are then obtained from the parameterization ( p3| ) and improved by 
further simulations, aiming at a precision of Ag^ ^0.1. This will be finished soon. 

Small mismatches of the simulation results w.r.t. the target values, i.e. g^{L\) = 4.484 and 
L\m\{L\) = 0, may be corrected by the non-perturbative jS-function and the mass derivative of the 



coupling [E6, 29] 



2.3.3 Parameters for HQET simulations in L| x T with T = L2,L2/2 

To prepare for the power-divergent subtractions in the volume of extent L2 = 2L\ 1 fm within 
the effective theory that eventually provide the link to HQET observables in the physically large 
volume (of extent Leo), the two-flavour theory will have to be simulated at typical resolutions of 
about L2/a = 8, 12, 16 and lattice spacings corresponding to 5.3 ^ j8 < 5.9. 

For fixing the necessary simulation parameters by means of the condition of fixed coupling 
g^{L2), one can rely again on the scale L* and its ratio to L2, r2 = L2/L*| continuum ~ 1-6, and infer 



the wanted pairs (L2/<3,j8) from eq. (2.5) 
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3. Computation of the missing improvement and renormalization factors 

We now present our non-perturbative determination of the improvement coefficient and the 
renormalization constant Z in the j3 -range relevant for the matching of HQET to QCD in small 
volume, such that the RGI heavy quark mass can be set to desired values z = L\M. 

Our generation of unquenched gauge configurations with SF boundary conditions for A^t = 2 



0(a) improved massless Wilson quarks employs the hybrid Monte Carlo (HMC) algorithm [ |30| ] in 



its variant used in the study of autocorrelation times in ref. It comprises multi-time-scale inte- 



gration schemes [ ]32| , |33| ] with mass preconditioning [g4j-^ on top of even-odd preconditioning. 



3.1 Non-degenerate current quark masses and estimators for b\ — bp, b^ and Z 

We proceed following the idea of imposing improvement conditions at constant physics, which 



was first advocated in [ ]27| ] and already applied to the present situation but for A^f = in []21|]. 

Since the definition (2.1) of Li via the renormalized coupling g^{Li/2) = g^{Lo), respectively 



the bare parameters in columns 2 - 4 of table |I| complying with it, have such a constant physics 
condition built in from the start, we can directly work at those pairs of (Lo/a,j3). With this as our 
choice of improvement condition, supplemented by the SF-specific settings of zero boundary fields, 
6 = 0.5 and — just for the purpose of this section — T /Lq = 3/2, the improvement coefficients 
bA — bp and b^ and the renormaUzation constant Z become smooth functions of gQ in the region 
where they are needed.^ 

Taking over any unexplained notations and details from refs. [^, 27 1 (and references therein). 



bf^ — bp, bm and Z can be determined by studying QCD with non-degenerate valence quarks. Treat- 
ing the latter in the quenched approximation, the structure of the 0(a) improved theory in conjunc- 
tion with a massless renormalization scheme retains the relative simplicity of the Nf = case elabo- 



rated in ref. [27|. For instance, the improvement of the off-diagonal bilinear fields X* =X^ ± /X^, 
X =A^,P, emerging as a consequence of the broken isospin symmetry in flavour space, is the same 
as in the degenerate case, except that the Z^-coefficients now multiply the average ^(awiq + amq,y) 
of the subtracted bare quark masses, m^j = mo,,- — nic, which themselves are separately improved 
for each quark flavour: 

mqj = niqj {\+bmamc^.i) . (3.1) 

(Here and below the indices i,j label the different quark flavours.) Identifying the valence flavours 
in the isospin doublet with a light (strange) and a heavy (bottom) quark, the corresponding PCAC 
relation reads 

d^A^ {x) = {mi + mj)P^ {x) , (3.2) 

and the renormalization constants Za and Zp that come into play upon renormalization are just 
those known in the theory with two mass-degenerate quarks. 

Accordingly, the SF correlation functions involving the axial current and the pseudoscalar 
density generalize to /'/{xq) = — ^(Aq (x)O^) and fp{xo) = — 5 (P+(x)0^), with pseudoscalar 
boundary sources decomposed as = ±iO^ where O" = a^Ly.z C(y)75 5 C(^)- Then the 



^Although the difference of coefficients bA — bp is actually not needed for fixing the RGI mass through eq. ( 2.2 1 
we include it in the present discussion. 
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improved bare PCAC (current) quark masses^ as functions of the timeslice location xq are given by 

mij{xo;LQ/a, T/Lq, 6) = jj— , (3.3) 

2fp'{xo) 

where only here we explicitly indicate their additional dependence on Lo/a, T /Lq and 6. In the 
degenerate case, / = j, the correlators assume the standard form, and m,-,- just reduces to the current 
quark mass of a single quark flavour that is prepared by a corresponding choice of equal values for 
the associated hopping parameters, k", = Kj. Also the precise definition of the lattice derivatives in 
eq. ( |33| ) matters. As it is written there, do = ^ {d^ + ^o) denotes the average of the ordinary forward 



and backward derivatives, but as in refs. [gjj, |27|] we have used the improved derivatives 

5o ^ 3o (1 - i a^d^o) , dtdo ^dtdo{l-^ a^dtdo) (3.4) 

as well, which (when acting on smooth functions) have 0(^Qa^,a^) errors only. 

For their numerical calculation, the coefficients bA — bp, bm and the finite factor Z = Z^Zp/Za 
(see eq. ( |2.3[ )) are isolated by virtue of the identity 



2 V'"q.' 



+ niqj) +jbm {am^j + am^j) - ^ (^a - ^p) a {mqj + mqj-) + O (a^) . (3.5) 



It is obtained by equating the expression for the 0(a) improved renormalized quark mass in terms 
of the bare PCAC mass with the alternative expression in terms of the subtracted bare quark mass. 
Forming ratios of suitable combinations of degenerate and non-degenerate current quark masses in 



the representation ( p.5| ) then enables to derive direct estimators for bx — bp, bm and Z [|27[]: 



Rap = 7 r-, r = 6a - op + O ( amq,i + amq,2 ) , (3.6) 

(mil -m22)(fl'mqj -fl'mq.2) 

4(mi2-m33) u i r,f , A n\ 

Rm = 7 77 T = byn + O [aniq i +amq 2) , (3.7) 

(mil -m22)(amqj -fl'mq,2) 

with mo,3 = ^(mo,i +mo^2)> neglecting other quark mass independent lattice artifacts of 0(a). For 
the renormalization constant Z an analogous formula holds even up to O(a^) corrections, 

Rz = — '^^^ + (bA-bp-bm)(amii+am22) = Z + O (a^) , (3.8) 

mqj-mq,2 

if the correct value for Z^a — ^p — = ^ap — (only involving correlation functions with mass 
degenerate quarks) is inserted. Note that generically the combination Z = Z^Zp/Za is a function 
of the improved bare coupling, = g^^l +bgamq). Since, however, we only consider light sea 
quarks that are massless (i.e. such that mi k, 0) and the valence quarks are anyway treated in the 
quenched approximation, this fact can be ignored here. 

Still, to complete our definition of the line of constant physics, values for the bare PCAC 



masses of the valence quarks must be selected. As in [|21[], we consider two pairs, 

choice 1 : LomnssO, Lom22~0.5, (3.9) 
choice 2 : LomnssO, Lom22~2.4. (3.10) 



■'This expression for the PCAC masses is only 0(a) improved up to a factor 1 + 5 (^a — bp) {aniq^i + am^j) for quark 
mass dependent cutoff effects. 
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set 


Lo/a 








Lom22 


bh-bp 




Z 


1 


10 


6.1906 


0.136016 


0.134318 


0.4929(6) 


-0.0006(9) 


-0.6643(9) 


1.1046(1) 




12 


6.3158 


0.135793 


0.134378 


0.4952(9) 


-0.003(2) 


-0.668(2) 


1.1050(2) 




16 


6.5113 


0.135441 


0.134387 


0.492(1) 


-0.006(2) 


-0.667(3) 


1.1044(2) 




20 


6.6380 


0.135163 


0.134356 


0.5005(9) 


-0.005(3) 


-0.669(3) 


1.1038(2) 


2 


10 


6.1906 


0.136016 


0.127622 


2.2909(5) 


+0.0727(4) 


-0.5655(3) 


1.0954(1) 




12 


6.3158 


0.135793 


0.128755 


2.3475(7) 


+0.0513(5) 


-0.5785(5) 


1.0974(1) 




16 


6.5113 


0.135441 


0.130146 


2.407(1) 


+0.0297(7) 


-0.5964(8) 


1.0995(1) 




20 


6.6380 


0.135163 


0.130965 


2.4433(8) 


+0.0215(6) 


-0.6076(8) 


1.1002(1) 



Table 1: Lattice parameters and numerical results on the improvement coefficients bp^ — bp and bm and on 
the renormalization constant Z. The parameters {Lo/a,l5 ,Ki) referring to the light (sea) quark sector have 
fixed SF coupling, g^{Lo) =2.9(1), and vanishing quark mass such as to meet the constant physics condition 
of Section 2. 1 . Our results for /ja — bp, b^ and Z are based on statistics varying from 0(300) measurements 
(Lo/a — 20) to 0(2000) measurements {Lo/a — 10). The upper set refers to "choice 1", eq. (3.9), where 
the heavy quark mass is kept at Lo'W22 ~ 0.5, while the lower set belongs to "choice 2" with Lo'«22 ~ 2.4, 
eq. (3.10). The condition Lquih « is fulfilled up to negligible deviations of about 0.015 at most. 



The first choice on Loni22 is motivated by the quenched investigation [|27p, where it was argued to 
be advantageous w.r.t. the size of 0(a) ambiguities encountered, while with the second choice one 
is closer to the typical b-quark region itself. Satisfying these conditions on L0/M22 for all (Lo/a, j8) 
in table |l] demands to properly adjust the hopping parameter, called Kj, above, that is responsible 
for the mass value of the heavy valence quark flavour. This in turn amounts to prior evaluations of 
the relevant correlation functions on the dynamical gauge background for some trial guesses of Kh, 
in order to estimate the heavy flavour's PCAC mass through eq. (|3.3[) and to tune it to the values 



dictated by eqs. (3.9) and ( p.lQ| ) up to a few percent.^ The resulting hopping parameters are given 



in the fifth column of table 
3.2 Results 



The technical aspects of the analysis to compute the estimators ( [3.6| ) - ( [3.8| ) from the numerical 
data on the heavy-light SF correlation functions by means of the PCAC masses niij for the various 
(degenerate and non-degenerate) valence quark mass combinations are the same as in refs. 
The correlators have been evaluated on our dynamical gauge field configurations, which were gen- 
erated on Lq X T lattices with T = 3Lq/2 and massless sea quarks (thus complying with the above 
requirement Kj = resp. Lomn ^ for the light valence quark flavour) and which were separated 
by 5 - 10 HMC trajectories of length one. As for the niij themselves, they have been calculated 



from the local masses, eq. ( |3.3[ ), using improved derivatives ( [3 .41 ) throughout and averaging over 
the central timeslices Lo/(2a),... ,{T — Lo/2)/a to increase statistics. Being secondary quantities 
in particular, the statistical errors of the masses and of the Rx, X = AP,m,Z, obtained from them 
were estimated by the F-method [B^, which directly analyzes autocorrelation functions. 



'^Similar to the situation in refs. |gjj,|27|], this is to sufficient precision equivalent to keeping fixed the corresponding 
renormalized masses LoZ^m/Zp, as for the considered couplings the entering renormalization constant barely varies. 
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Figure 2: Two sets of non-perturbative results for bA — bp in the considered region of bare couplings, 
referring to our two choices of quark masses, together with the prediction from one-loop perturbation theory. 
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Figure 3: The same as in figure ||but for the improvement coefficient bn 



Our non-perturbative results on bp 



bp, bm and Z are also listed in table [T[ As a consequence 

AP,m,Z, become smooth 



of the underlying constant physics condition ( [2.1[ ), the estimates Rx, X 
functions of the bare coupling, = 6/j3 . This is well reflected in figures ^ - Q where our results 
are shown in comparison with the one-loop perturbative predictions 

The overall ^Q-dependence of our results is qualitatively similar to the quenched study [ 21 ] 



and even comparable on the quantitative level. Whereas Rap is compatible with a nearly vanishing 
bA — bp, as predicted by leading-order perturbation theory, for "choice 1" of quark masses and 
appears to approach this Une quite rapidly as — for "choice 2", one observes for both choices 
significant deviations of the sets of estimates for b^ and Z from the leading perturbative behaviour 
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Figure 4: The same as in figure g but for the renormalization constant Z. 
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Figure 5: Left: Difference of the two sets of results in table |l] on the renormalization constant Z versus 
(fl/Lo)^- Right: The same for the improvement coefficient bj^ where, however, the ambiguity inherent in any 
improvement condition imposed is of 0(fl). The open black triangles display the corresponding quenched 
results from [ pl] | for comparison. 



in the weak coupling region considered. Since one expects the perturbative curves eventually to be 
approached in the limit — > also in case of Z>m and Z, the curvature seen in our numbers hints 
at a more complicated structure of (unknown) higher-order terms. Hence, we have to conclude 
that if an improvement condition were used in a region of stronger couplings, which would no 
doubt lead to a rather different set of data points, simple one-loop perturbation theory would not 
be an adequate guide for the continuation of and Z to weak couplings. On the contrary, this 
would induce a source of uncertainty in results deriving from them that is difficult to control and, 
therefore, highlights the importance of employing improvement conditions in the jS-range relevant 
to the actual apphcation. 

Of course, any other estimate Rx (i.e. stemming from a different choice of renormaliza- 
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tion/improvement condition) may yield a different functional dependence upon gg, but its differ- 
ences are again smooth functions that vanish in the continuum limit with a rate proportional to a/Lo 
(for improvement coefficients) or {a/Lof- (for renormalization constants). These intrinsic 0(fl'") 
ambiguities (« = 1 , 2) imply that rather than a numerical value at some given j3 , the important infor- 
mation lies in the correct gQ-dependence of the estimators Rx, X = AP,m,Z, obtained at constant 
physics. To demonstrate this, we also investigated a few alternative improvement conditions, which 
are either provided by defining the estimators Rx with standard instead of improved derivatives or 
by the two quark mass choices, eqs. ( |3.9[ ) and ( p.lO| ), themselves. As an example we plot in the 
left panel of figure |5| the difference AZ{gl) = Z(gg)| choice i -^(gQ)! choke 2 versus (a/Lo)^, which 
clearly exhibits a linear approach towards zero. Other cases behave similarly, e.g. the 0{a) ambi- 
guities for Abmigl) = (go) I choice 1 " feo ) I choice 2 in the right panel of figure | are found to be 
quite small, and their magnitude rapidly decreases as tz/Lo ^ 0. 



4. Outlook 

Apart from the elements sketched at the end of Section |2[ which partly are already in progress, 
the computation of the b-quark mass at the 1/m-order of HQET along our strategy illustrated in 
figure |I] still requires N{ = 2 simulations in L2 ~ 1 fm as well as in physically large volumes of 
about Loo > 2fm. Particularly for the latter we plan to switch to QCD with periodic boundary 
conditions and to use the technique of low-mode deflation [|^] in connection with all-to-all quark 



propagators [ |40| ] for the numerical evaluation of correlation functions. 



As further interesting directions for future work let us mention the non-perturbative tests of the 



HQET expansion in the spirit of ref . 1 22 ] and the extension of our determination of improvement 



coefficients and Z-factors to the parameter range relevant for (large volume) charm physics. 
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